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Abstract. It is well appreciated that the description of overadiabatic convection affects the structure of the envelopes of lumi- 
nous asymptotic giant branch (AGB) stars in the phase of "hot bottom burning " (HBB). We stress that this important uncertainty 
in the modeling plays a role which is much more dramatic than the role which can be ascribed, e.g., to the uncertainty in the 
nuclear cross-sections. Due to the role tentatively attributed today to the HBB nucleosynthesis as the site of self-enrichment 
of Globular Clusters stars, it is necessary to explore the difference in nucleosynthesis obtained by different prescriptions for 
convection. We present results of detailed evolutionary calculations of the evolution of stars of intermediate mass during the 
AGB phase for the metallicity typical of the Globular Clusters that show the largest spread in CNO abundances (Z ~ 10~ 3 ). 
We follow carefully the nucleosynthesis at the base of the external convective region, showing that very different results can 
be obtained according to the presciption adopted to find out the temperature gradient within the instability regions. We discuss 
the uncertainties in the yields of the various chemical species and the role which these sources can play as polluters of the 
interstellar medium. 
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1. Introduction 

The intermediate mass stars (i.e. stars with initial masses 1 < 



M/M G 



7, hereinafter IMS) experience a phase of ther- 



mal pulses (TPs) shortly after the exhaustion of central he- 
lium (Schwarzschild & Harm 1965, 1967; Iben 1975, 1976). 
A CNO-burning shell supplies the global nuclear energy for 
~ 90% of the time; periodically, a thermal instability (asso- 
ciated with the narrow dimensions of the helium-rich layer 
above the CO degenerate core and with the steep dependence 
of the cross-section of 3c reactions on the temperature) triggers 
a runaway-like ignition of helium burning, which temporarily 
extinguishes the H-burning shell. As the mass of the CO core 
increases, the CNO-burning shell becomes hotter and hotter, 
so that the luminosity of the star increases: on the HR dia- 
gram the evolutionary tracks climb the so called "Asymptotic 
Giant Branch" (hence the term "AGB evolution", commonly 
adopted). During their evolution these structures suffer strong 
mass losses, which eventually stop their AGB evolution, leav- 
ing just a central remnant which later evolves into a CO white 
dwarf. 

In the last decades the interest in the detailed evolution of 
these stars has grown, also because they have been suggested 
as pollutors of the medium during the early phases of evolution 
of globular clusters (GCs), and responsible for the chemical 
anomalies (in terms, e.g., of oxygen-sodium and magnesium- 
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aluminium anticorrelations) which are observed in giants and 
turn-off stars of GCs. At least in the latter because of the low 
temperatures involved the chemical anomalies cannot be due to 
"in-situ" processes. The idea behind this hypothesis (to which 
we usually refer as "self-enrichment") is that the base of the 
external envelope of the most massive IMS during the AGB 
evolution may become so hot (7b ce > 30 ■ 10 6 K) that it triggers 
strong nucleosynthesis; the consequent changes in the origi- 
nary chemical abundances might easily reach the surface of the 
star because of the rapidity of the convective motions, so that 
the ejecta of these stars are contaminated by such an advanced 
nucleosynthesis. The velocity of the ejected material should be 
sufficiently low to survive ejection from the cluster (D'Antona 
et al. 1983; Ventura et al. 2001, 2002). 

The self-pollution scenario is not without its shortcomings, 
as pointed out by Denissenkov & Herwig (2003): the temper- 
atures needed to deplete oxygen would also lead to sodium 
destruction, and to magnesium isotopic ratios 25 Mg/ 24 Mg and 
26 Mg/ 24 Mg much larger than observed. 

Unfortunately, the theoretical models for the AGB evolu- 
tion are characterized by several uncertainties associated with 
the input physics which is commonly adopted in the stellar evo- 
lutionary codes: in the last decade the attention of modelists has 
been mainly focused on: a) the amount of extra-mixing from 
the base of the external region (Herwig 2000, 2004), which 
favours a larger efficiency of the so called third dredge-up phe- 
nomenon (i.e. the external convective region which reaches 
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stellar layers previously touched by efficient 3a burning after a 
thermal pulse); b) the rate of mass loss suffered by these stars; 
c) the nuclear network adopted, particularly concerning the nu- 
cleosynthesis of Neon, Sodium and Magnesium isotopes. 

The present work does not discuss the plausibility of the 
self-enrichment scenario, but is aimed at pointing out quantita- 
tively that, before going into the details of the chemical com- 
position of the ejecta of AGBs, it is mandatory to understand 
to which extent the results obtained depend on the assumptions 
made in calculating the models. 

We focus our attention on one of the most relevant uncer- 
tainties connected with stellar evolution, i.e. the treatment of 
convection, whose influence upon AGB evolution in the con- 
text of HBB is well documented in the literature, starting with 
the classic Renzini & Voli (1981) paper, and with the discussion 
by Sackmann & Boothroyd (1991) and Blocker & Schonberner 
(1991). Finally, D'Antona & Mazzitelli (1996) discussed that 
the "Full Spectrum of Turbulence" (FST, Canuto & Mazzitelli 
1991) model for convection leads to a much more efficient 
HBB than the Mixing Length Theory (MLT). Neverthless, in 
these last years, no mention has been made in many recent de- 
tailed computations of the description adopted for convection, 
and, more importantly, of the resulting uncertainty in the nu- 
cleosynthesis. 

We will compare the results obtained with the two local 
models currently available, i.e. MLT and the FST convective 
model; within the MLT framework, we discuss the role played 
by the free parameter a. 

We show that the main physical properties of the models 
depend strongly on the efficiency of convection, and discuss 
the implications for the self-enrichment scenario. 

2. The evolutionary code 

The stellar evolution models discussed in this paper were cal- 
culated by the code ATON, a full description of which can be 
found in Ventura et al. (1998) (ATON2.0 version). In this pa- 
per the interested reader may find a detailed description of the 
numerical structure of the code, and of the macro- and micro- 
physics used to simulate the stellar evolutions. 

The code has now been updated concerning the nuclear 
network, which has been widened in order to include more 
chemical elements and nuclear reactions, according to the de- 
tailed description given below. The current version is therefore 
ATON2.1. 

Here we briefly recall the most important input data 
adopted. 

2.1. The convective model 

The code allows us to calculate the temperature gradient within 
instability regions either by adopting the traditional MLT 
(Vitense 1953; Bohm-Vitense 1958), or the FST model (Canuto 
& Mazzitelli 1991; Canuto et al. 1996). The interested reader 
may find a detailed description of the physical differences be- 
tween the two models in Canuto & Mazzitelli (1991). 

Briefly, we recall that within the MLT scheme both the di- 
mensions of the convective eddies and the mixing length are 



assumed to be directly proportional to the local value of H p , 
the pressure scale height (I = aH p ); the free parameter a is 
calculated in order to reproduce the evolution of the Sun. The 
most recent estimates give a = 1.7. In the FST model the mean 
dimension of the convective eddies is found by integrating over 
the whole spectrum of the dimensions, and the mixing length 
is simply assumed to be I = z, the distance from the closest 
convective boundary. 

The differences between the results provided by the two 
models can be summarized as follows: 

- within high-efficiency convective regions (i.e. stellar re- 
gions where most of the energy is carried by convection) 
the level of overadiabaticity V - V a( j, (where V = 
required by the FST models is lower 

- within low-efficiency instability regions (e.g. low density 
convective zones) the FST model exhibits a superadiabatic- 
ity peak which is not found within the MLT framework 

2.2. Nuclear network 

With respect to Ventura et al. (1998) the nuclear network 
has been widened and now includes 30 chemicals:rm//,D, 
3 He, 4 He, 7 Be, 7 Li, 12 C, 13 C, 14 N, 15 N, 16 0, 17 0, 18 O, 18 F, 19 F, 20 Ne, 
21 Ne, 22 Ne, 22 Na, 23 Na, 24 Mg, 25 Mg, 26 Mg, 26 Al, 27 Al, 28 Si, 29 Si, 
30 Si, 31 P,n. The nuclear reactions considered are 64. They in- 
clude all the main reactions of the p-p, CNO, Ne-Na and Mg- 
Al chains, and the a captures of all nuclei up to 26 Mg. The 
relevant cross-sections can be taken either from Caughlan & 
Fowler (1988) or from Angulo et al (1999). 

2.3. Convective mixing and nuclear burning 

Since we expect a non-negligible fraction of the global nuclear 
release to be generated within the convectively unstable exter- 
nal regions, we decided to adopt for all the evolution models 
presented here a diffusive algorithm to deal with nuclear burn- 
ing within convective regions, in which nuclear burning and 
mixing of chemicals are coupled self-consistently. We there- 
fore solve for each element the diffusion equation (Cloutman 
&Eolll976): 

(f) = (f) nucl ^K 2 p)^S] a> 

stating mass conservation of element ;. The diffusion coeffi- 
cient D is taken as 

D = \ul (2) 

where u is the convective velocity and I is the convective scale 
length. 

Within this diffusive framework it is necessary to spec- 
ify the way with which convective velocities decay outside 
the convective boundaries (Deng at al. 1996a,b; Herwig et al. 
1997; Ventura et al. 1998). In agreement with Xiong (1985) and 
Grossman (1996) and supported by the numerical simulations 
by Freytag et al. (1996), we assume that convective velocities 
decay exponentially outside the formal convective boundary as: 

u = u » ex P ± ^ 
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Fig. 1. The tracks in the HR diagram of three models with ini- 
tial mass 5M calculated with three different tretaments of con- 
vection. Full: FST models; dotted: MLT model with a — 1.7; 
dashed: MLT model with a = 2.1 

where ut, and Pt, are, respectively, turbulent velocity and pres- 
sure at the convective boundary, P is the local pressure, £ a free 
parameter connected with the e-folding distance of the decay, 
and /thick is the thickness of the convective regions in fractions 
of H p . 

3. Early evolutionary phases 

We computed three evolutions starting with initial masses M = 
5M Q . With the exception of the convective model, all the phys- 
ical and chemical input data adopted are the same, i.e.: 

- Since our purpose is to test how the uncertainties connected 
with convection may reflect on the plausibility of the "self- 
enrichment" scenario, we adopt a chemistry which is typi- 
cal of those GCs which show most of self-enrichment like 
M13, NGC6752 (Gratton et al. 2001, Sneden et al. 2004), 
i.e. Z = 0.001 and Y = 0.24. For all chemicals other than 
helium we adopted solar-scaled abundances. 

- We assumed extra mixing from the external border of any 
convective region as given by a parameter ( - 0.02 for 
the exponential decay of velocities starting from the for- 
mal convective border; this is in agreement with the cali- 
bration given in Ventura et al. (1998). Since our purpose is 
to discuss the uncertainties connected with the treatment of 
convection, no inwards overshooting from the base of the 
external convective region was adopted. 

- For mass loss, Blocker's (1995) prescription was used in 
all cases: the free parameter 77s entering Blocker's formula 
was set to tjr = 0.02, in agreement with the calibration for 
mass loss during the AGB evolution given in Ventura et al. 
(2000). 

- The Angulo et al. (1999) cross-sections were used. 

- Two models were calculated with an MLT treatment of 
convection with parameters a — 1.7 (MLT17 model) and 
a = 2.1 (MLT21 model). We also calculated an evolution 
model with the FST convective prescription (FST model). 

Fig- HI shows the tracks in the HR diagram corresponding 
to the three evolution models discussed above. We note that the 
tracks are very similar during the two major phases of nuclear 
burning within the central regions, because convection within 
burning cores is so efficient that the gradient is practically adi- 
abatic, independently of the convective model. The only differ- 
ences can be noticed in the colors of the red giant branch (RGB) 
and the early AGB evolution. These are the regions of the HR 
diagram where the stars develop a very extended convective en- 
velope, in which a great part of the energy is carried by radia- 
tion because of the low convective efficiency: in these cases the 
temperature gradient from the CNO-burning shell to the sur- 
face (hence the effective temperature of the star) depends sen- 
sibly on the convective model adopted. A full discussion of the 



Fig. 2. The variation with time of the temperature at the base of 
the convective envelope (top), luminosity due to helium burn- 
ing (middle), and core mass (bottom) of a 5M FST model dur- 
ing the AGB evolution. Time was counted from the beginning 
of the AGB evolution. 

Fig. 3. The variation with time of luminosity (top), mass loss 
rate (middle) and total mass (bottom) of the model presented in 
Fig.E 

Fig. 4. Evolution with time of surface 3 He (top), lithium (mid- 
dle) and fraction of nuclear energy due to proton-proton reac- 
tions (bottom) of the FST model presented in Fig.|2] 

Fig. 5. Evolution with time of the surface chemical abundances 
of the CNO elements (top) and of some neon, sodium and mag- 
nesium isotopes (bottom) of the FST model presented in Fig. 13 

dependency on the convective model of the effective tempera- 
ture of IMS during various evolutionary phases can be found in 
Ventura & Castellani (2004). 

What is most interesting here is that the main physical prop- 
erties and the duration of the various evolutionary phases are 
approximately the same for the three models, because in any 
case the nuclear sources are either well within regions of ra- 
diative stability, or inside central cores where the convective 
efficiency is extremely high, making the gradient almost adia- 
batic. 

The phase of H-burning lasts 90Myr; during it the star 
develops a central convective core of ~ 1.4M , which pro- 
gressively shrinks in mass until H-burning is extinguished in 
the central regions. Soon after hydrogen exhaustion the star 
expands, and surface convection reaches inner layers which 
were previously at least partly touched by nuclear burning (first 
dredge-up). The lowest point (in mass) reached by the base of 
the external envelope is for all the three models ~ 2.1M away 
from the centre; the whole duration of this phase is ~ 70, 000 
yr. 

Shortly afterwards the central core becomes hot enough to 
ignite 3a reactions. During helium burning the star develops a 
central convective core of ~ O.5M ; the total duration of this 
phase for the three models discussed is ~ 12Myr. 

After the central helium exhaustion, 30- reactions operate 
in an intermediate layer above the core of carbon and oxygen. 
All the stellar layers above this He-burning shell expand, and 
the CNO-burning shell is temporarily extinguished. Convection 
penetrates inwards, reaching a layer which is O.9M away from 
the centre (second dredge-up). The total duration of this phase 
of inner penetration of surface convection is 0.55 Myr, after 
which the CNO-burning shell is activated again, and the star 
begins the AGB evolution. 

4. Standard AGB evolution 

Shortly after the second dredge-up the external layers of the 
star contract and heat, thus favouring the re-ignition of the 
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CNO-burning shell. Hereinafter, the CNO cycle will be the 
only nuclear source supporting the star, with the only excep- 
tion of the contribution of a 3a burning shell which once every 
~ 3000 yr is ignited in a thermally unstable way, leading to a 
thermal pulse (TP). 

Fig. I3- 14] show the temporal evolution of some physical 
quantities related to the evolution of the FST model with initial 
mass 5M Q . We note a rapid increase of the stellar luminosity 
(Fig. top panel), which is associated with the fast increase 
of 7bce (Fig. top panel), which is almost doubled after three 
TPs. We will see that the rapidity of this rise is strongly depen- 
dent on the convective model adopted. 

The periodic strong drop of the global luminosity and of 
7bce marks the ignition of TPs, which trigger an expansion of 
the layers above the helium burning shell with the consequent 
stop of CNO burning and shrinking of the surface convective 
envelope. 

The CNO luminosity is generated in a shell which is 6M ~ 
1O~ 4 M wide, and which progressively moves outwards, burn- 
ing at higher and higher gravities and temperatures: this de- 
termines the increase of luminosity shown in the upper panel 
of Fig. [5] which is eventually halted by the strong mass loss, 
which reduces the mass of the envelope, as can be seen in 
the bottom panel of Fig. [3] In the same figure we see the 
rapid increase of the mass loss with luminosity (middle panel): 
within Blocker (1995) prescription, and with the adopted value 
of the free parameter t]r, we see that a maximum value of 
~ 2 x 10~ 4 M o /yr is attained, almost in conjunction with the 
maximum luminosity of the star. Therefore, the rapidity with 
which the luminosity increases plays a fundamental role in 
determining the total duration of the whole AGB evolution, 
triggering a fast reduction of the mass of the envelope. For 
the model discussed here, we see that the AGB phase lasts 
~70,000yr. 

For a few TPs we note from the bottom panel of Fig.|4]that 
a non-negligible fraction of the nuclear energy release is due 
to the proton-proton (p - p) chain: during the 4th interpulse 
phase this latter channel provides ~ 10% of the luminosity of 
the star. This is due essentially to 3 He burning at the base of the 
external envelope (see the top panel of Fig.|4}, which becomes 
efficient as soon as the temperature reaches ~ 4 x 10 7 K. The 
p - p contribution becomes negligible as soon as the surface 
3 He is all burnt (bottom panel). We also note that this is asso- 
ciated with the lithium production via the Cameron & Fowler 
(1971) mechanism, as can be seen by the large increase of the 
surface lithium abundance shown in the middle panel of the 
same figure. The end of 3 He burning at the base of the convec- 
tive envelope is associated with the temporary plateau in the 
global luminosity rising which can be detected in the top panel 
of Fig. |3] after - 12000 yr. 

We end this general description of AGB evolution with a 
glance at the nucleosynthesis at the base of the convective en- 
velope: we recall that convection is generally extremely fast in 
homogenizing the whole external region; therefore the ejecta 
of these stars are practically determined by the changes of the 
chemical abundances at the base of the external zone. 

We note from the top panel of Fig. [5] that for the first 
~ 15000 yr the oxygen abundance is unchanged, while 14 N 



Fig. 6. The comparison of the temporal variation of the lumi- 
nosity of three models with initial mass M = 5M Q calculated 
with three different prescriptions for the treatment of the re- 
gions unstable to convection. 

is produced at the expenses of I2 C: this is the signature of a 
CN cycle, which is operating at temperatures not exceeding 
~ 7 x 10 7 K. Later on, when 7b ce exceeds ~ 8 x 10 7 K, oxy- 
gen is depleted at the base of the envelope, with the consequent 
increase of the I4 N and 12 C abundances. After ~ 15 TPs, we 
can cleary see the signature of the Third Dredge-Up (TDU) 
in the strong increase of the surface 12 C abundance following 
each TP. More precisely, we find the first TDU episode after 
13 TPs; the efficiency parameter A (defined as the ratio be- 
tween the mass which is dredged up following each TP and 
the mass gained by the H-exhausted region since the precedent 
TP) reaches an asymptotic value of A ~ 0.5 after 19 TPs. 

These results confirm, at least within the FST framework, 
the possibility of attaining deep oxygen burning at the base of 
the external convective zone of these stars: the final I6 abun- 
dance is > 10 times lower, while the ejecta would have an 
average oxygen content a factor of ~ 3 lower than the initial 
abundance. 

From the bottom panel of Fig. [5] we see an early phase 
of 22 Ne destruction, which favours a temporary production of 
sodium, increased by another factor of ~ 2 (in the comparison 
of the current sodium abundance with the initial value, we must 
recall that the second dredge-up had already favoured a sodium 
increase of a factor ~ 3). Shortly after the beginning of oxygen 
depletion, the surface sodium abundance starts decreasing due 
to the NeNa chain. 

Later on, when TDU becomes operative, a considerable 
amount of 22 Ne is carried outwards to the surface following 
each TP: this favours sodium production again, as we see in 
the bottom panel of fig|5] The overall sodium abundance of the 
ejecta would be sligtly lower than the initial value. 

Finally, we note that 24 Mg is also depleted at the base of 
the external convective zone: the final abundance is lower with 
respect to the initial value by a factor of ~ 20. 

5. The influence of the convective model 

We stressed in Sect. 3 that prior to the TPs phase the three 
evolutions are very similar. Here we focus our attention on the 
AGB phase, keeping in mind that at the beginning of this phase 
the structures of the three models are practically identical. 

In Fig. [6] we show the variation with time of the luminos- 
ity for the three models discussed. We see a large difference 
starting already at the first TPs; the more efficient the convec- 
tive model is, the more rapid is the increase of luminosity as the 
evolution proceeds. We see that the FST model achieves a max- 
imum luminosity which is ~ 70% larger than the corresponding 
value of the MLT17 model. 

This behaviour can be seen as an extreme case of the 
break-down of the core mass vs. luminosity relation which 
occurs when the hydrogen burning region is not detached 
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Fig. 7. The variation with time of the total mass of the models 
shown in Fig.|6] 

Fig. 8. The variation with time of the total luminosity (top) and 
of the maximum CNO nuclear energy release (bottom) for the 
three models discussed in Fig. [6] Only the early evolution along 
the AGB is shown. 

Fig. 9. The internal nuclear structure of the models presented in 
Fig-Elat the maximum luminosity between the 3rd and 4th TPs. 
The abscissa gives the distance (in solar masses) from the cen- 
tre of the star, and the ordinate shows the coefficient for nuclear 
energy release. The thin vertical lines indicate the location of 
the inner border of the convective envelope. 

from the convective region (Tuchman et al. 1983, Blocker & 
Schonberner 1991). 

Such a large difference in the evolution of luminosity is in- 
evitably reflected in the mass loss rate: we see from Fig. a 
much faster consumption for the FST model, which loses the 
whole envelope within ~ 70,000 yr, to be compared to the 
~ 130, 000 yr associated to the MLT21 model, and ~ 200, 000 
yr of the MLT17 model. We therefore can see that the duration 
of the whole AGB phase may differ by a factor of ~ 3 according 
to the convective model adopted. 

Before entering into the details of the consequences which 
these differences may have for the chemistry of the ejecta, 
hence for the way in which these sources pollute the interstellar 
medium, we examine the reasons for such large differences. 

We show in Fig. [8] the comparison between the evolution 
of the three models, limited to the first 10,000 yr of AGB, or, 
equivalently, the first 3-4 TPs. The upper panel shows the total 
luminosity of the star, in the bottom panel we report the maxi- 
mum value of the CNO-burning coefficient for nuclear energy 
release. While the latter quantity is very similar for each of the 
three models, the luminosities diverge. We may therefore con- 
clude that while the internal structure of the models in terms of 
thermal stratification and chemical profile is similar up to the 
CNO-burning shell, some difference must arise in the immedi- 
ate proximity of the latter, close to the inner boundary of the 
convective envelope. 

To better clarify this point, we show in Fig. |9] the interior 
stratifications of the coefficient for the nuclear energy release; 
for all the three models the figure refers to the point of maxi- 
mum luminosity during the fourth interpulse phase, ~ 9000 yr 
after the beginning of the AGB evolution. The mass coordinate 
of the MLT17 and MLT21 models has been artificially shifted 
in order to have the peak value of e nuc i at the same abscissa. 
This shift is much less than 1O~ 3 M . The thin vertical lines in- 
dicate the location of the bottom of the convective envelope for 
each model. 

For each of the three models we see a very similar be- 
haviour. In the proximity of the e nuc i peak there is a secondary 
maximum at the border of the surface convection, which carries 
3 He-rich material into the CNO-burning shell; during these first 
TPs, as already discussed in Sect.4, the p - p contribution is 



non-negligible (see the bottom panel of Fig.|4j. Also, we see a 
secondary maximum ~ 2 x 1O~ 4 M away, which is entirely due 
to lithium burning; the contribution of the latter to the overall 
energy release, however, is negligible. 

If we compare the three models, we note a strong similar- 
ity up to the peak of the shell and just beyond, with the only 
difference that the location of the inner border of the convec- 
tive external zone is closer to the CNO peak within the FST 
model, while it is a few 10~ 5 M o further away in the MLT17 
case; the MLT21 model is intermediate. The proximity of the 
convective border to the CNO peak carries more 3 He within a 
region whose temperature is typical of CNO-burning regions 
(T ~ 7 x 10 7 K), thus triggering a considerable extra luminos- 
ity. This extra contribution is smaller the more distant the CNO 
peak and the base of the envelope are, and it is practically neg- 
ligible within the MLT17 model. 

From the above discussion it becomes evident that in un- 
derstanding the different luminosities of the models, despite 
the similarity characterizing their structures up to the CNO- 
burning shell, a key role is played by the exact location of the 
inner border of the convective envelope, or, equivalently, its 
distance from the peak of the CNO shell. Within the frame- 
work of the Schwartzschild's criterion, which is used to fix the 
neutrality point according to the condition V ra d = V a d, the loca- 
tion of the border is furhter inward the steeper the V la d profile 
towards the external layers. 

The comparison of the variations of the V ra( j within the 
three models presented in Fig. |5]is reported in the left panel 
of Fig. [TO] We see in all cases that V ra( j increases from the peak 
of the shell, because of the rapid decrease of temperature and 
pressure and the increase of the opacity; near the convective re- 
gion the FST V ra d profile keeps the same slope, while the MLT 
models, particularly MLT17, become progressively flatter, thus 
moving the neutral point, where V ra d = V a d, outwards. 

The reason for this behaviour can be found in the difference 
in the efficiency of the three convective models, and the relative 
effect on the temperature profile. It is well known that the FST 
model is more efficient in high-efficiency convective regions, 
thus requiring a lower degree of overadiabaticity to carry the 
same amount of energy flux by convection: in the layers next 
to the inner border of surface convection shown in Fig. [9] the 
FST overadiabaticity is ~ 0.8 dex smaller. In a small region 
very close to the border the trend is reversed, because the MLT 
models attain an extremely low level of overadiabaticy: this is 
due to a physical inconsistency of the MLT assumption that the 
mixing length is proportional to H p , which is meaningless near 
the border, where the same quantity is expected to vanish, as is 
consistently described in the FST. 1 

Within the instability regions, the fraction of energy which 
is carried by radiation decreases more rapidly in the FST than in 
the MLT models; the necessity of keeping the radiative flux ,F ra d 
higher requires a lower and flatter profile of V ra d (we recall that 
^rad K 1 /V ra d; it is the necessity of matching this profile within 
the convective region which forces the MLT radiative gradient 

1 yet this region is so efficient convectively that this has no conse- 
quences for the thermodynamical structure of these regions, the gradi- 
ent being very close to the adiabatic value. 
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Fig. 10. Internal distribution of radiative gradient (left), and temperature (right) for the same models shown in fig|9] 

Fig. 11. The variation of the abundances of some chemicals at the surface of the star as a function of the total mass (decreasing 
due to mass loss). Left: CNO abundances; Right: Neon, Sodium and Magnesium abundances. 

profile to bend near the border, thus shifting the location of the Table 1. Chemistry of the ejecta. 



neutrality point outwards. 

The higher efficiency of the FST model thus has two im- 
portant consequences: 

- The FST V la( j profile is steeper, favouring a more internal 
location of the formal border; this in turn favours a larger 
luminosity, because we have an efficient extra contribution 
of the p — p chain. 

- The lower value of the overadiabaticity required in the 
FST model (and, partially, in the MLT21 model) leads to 
a smaller temperature gradient within the instability region 
(see the right panel of Fig. I10>. which, in turn, acts in favour 
of a larger contribution of the extra luminosity supplied by 
3 He burning within these stellar layers. 

From a chemical point of view the differences among the three 
models are also relevant. Fig. ^2 shows the variation of the 
abundances of some key elements included in our nuclear net- 
work; in the left panel we report the evolution of the CNO sur- 
face abundances, in the right 2() Ne, 23 Na and 24 Mg are shown. 
In this case, to have an idea of the mean chemical composition 
of the ejecta of the star, we use the mass as abscissa instead of 
time (Ventura et al. 2001); the difference is relevant, because 
most of the mass is lost in correspondence of the maximum 
luminosity, in a time interval which is very short. 

We note in the left panel of Fig. \^\ that all three models 
achieve oxygen depletion, though in MLT17 the final depletion 
is just by a factor of ~ 2. The MLT17 model, which yields a 
longer lifetime, experiences more TDUs, thus it produces more 
12 C and 14 N. As in the FST case, we find that the parameter 
connected to the efficiency of the TDU approaches an asymp- 
totic value of A ~ 0.5 after ~ 20 TPs. 

In the right panel we note the different behaviour of sodium. 
In all three cases we have an early phase of production (which 
is due to the combined effects of the second dredge-up and to 
the destruction of 22 Ne at the beginning of the AGB evolution), 
and a later phase of depletion; in the MLT models, however, 
the effect of several TDUs make the surface content of sodium 
rise again, and reach extraordinarily large values, particularly 
within the MLT 17 model. We also note a different degree of 
24 Mg depletion in the three cases, 24 Mg reduction being larger 
the larger the efficiency of the convective model. 

Table [^summarises the chemical composition of the ejecta 
for each of the three models discussed. As we may see there are 
just two robust conclusions which we can draw, independently 
of the convective model adopted: 

- The fraction (in mass) of helium in the expelled material 
is Y ~ 0.31 in all cases: this can be understood on the ba- 
sis of the fact that it occurs essentially during the second 
dredge-up, while the following production during AGB is 
negligible. 



Model 


Y{a) 


W(b) 


[C12](c) 


[M4](c) 


[0161(c) 


[Na23](c) 


[Afg24](c) 


MLT 11 


0.32 


2.0 


-0.3 


1.51 


-0.40 


0.78 


-0.45 


MLT21 


0.32 


1.9 


-0.5 


1.34 


-0.55 


0.60 


-0.83 


FST 


0.31 


1.9 


-0.7 


1.13 


-0.60 


-0.16 


-0.95 



a Helium mass fraction of the ejecta. 

b Lithium content of the ejecta on the scale log(/V( 7 Li)) = 

log( 7 Li/H) + 12.00. 
c Average abundances with respect to the initial value: [X] = 

log^ejecta/^mitial)- 



- The lithium content of the ejecta is log(e( 7 Li)) ~ 1.9, 
slightly less than the population II average abundance. This 
is because the phases of lithium production and destruction 
happen at the very beginning of the AGB evolution, when 
the differences among the three models are small. 

For all the other elements, the pollution of the interstellar 
medium is extremely dependent on the convective model: 

- In the MLT 17 case we would expect little oxygen deple- 
tion, and extremely large sodium and nitrogen enhance- 
ment. The sum C+N+O would be increased by ~ 0.8dex, 
due to the effects of several TDUs. We would also prob- 
ably expect strong s-process enrichment. The isotopic ra- 
tios of magnesium 25 Mg/ 24 Mg and 26 Mg/ 24 Mg would be 
slightly less than unity. These results are consistent with 
recent AGB models of the same metallicity presented by 
Fenner et al. (2004). 

- In the MLT2 1 model we would expect a larger oxygen de- 
pletion and lower sodium enhancement. The increase of the 
sum C+N+O would be lower (~ 0.4dex). The isotopic ra- 
tios of magnesium would both be around ~ 3. 

- The FST model achieves the largest oxygen depletion, 
while the sodium content is slightly underabundant with 
respect to the solar value. The C+N+O abundance is con- 
served; the isotopic ratios of magnesium are similar to the 
MLT21 case. 

In summary, we have seen that the stellar yields depend 
dramatically on the convection model. There are many other 
uncertainties in the results, but those are less critical. 

We take as an example the sodium content of the ejecta: we 
saw that the results change significantly according to the effi- 
ciency of convection at the base of the outer convective zone. 
Fig. ^] shows the variation with mass of the surface sodium 
abundance within the three models discussed so far (full, dot- 
ted, and dashed lines), plus a further FST model (dashed-dotted 
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line) calculated by assuming the lower limit given in Angulo et 
al. (1999) for the reactions destroying sodium ( 23 Na(p, y) 24 Mg 
and 23 Na(p, a-) 20 Ne) and the upper limit for reactions creating 
sodium ( 22 Na(p, y) 23 Mg and 22 Ne(/?,y) 23 Na). We see that the 
influence of changing all the cross-sections involving sodium 
is far less than the role played by convection. 

The problems of reaction rates and mass loss will be dis- 
cussed in a forthcoming paper. 

6. What implications are there for the 
self-enrichment scenario? 

Deep spectroscopic observations of giant and TO stars within 
globular clusters have shown star-to-star variations in surface 
chemical composition. Particularly, the observations of NGC 
6752 by Gratton et al. (2001) showed the presence of an O-Na 
anticorrelation in all stars observed, independently of their evo- 
lutionary stage, including those near the TO. Within the same 
cluster, strong hints of a Mg-Al anticorrelation were also found. 
Similar results were obtained for other clusters (NGC 6397, 
M30 and M55: Carretta 2003; M71: Ramirez & Cohen 2002; 
M5: Ramirez & Cohen 2003). These results have also been re- 
cently confirmed by a deep analysis of the surface abundances 
of giant stars belonging to M3 and M13 by Sneden et al. (2004), 
where a clear C-O anticorrelation was found. 

The aforementioned results indicate that the surface mat- 
ter of these stars has been subject to nuclear processing though 
CNO, Na-Ne and Mg-Al chains. Since these anomalies were 
observed also in TO stars, where the internal temperatures 
are too low to allow such nuclear reactions, it was suggested 
(Ventura et al. 2001; Ventura et al. 2002) that the stars dis- 
playing such surface chemical anomalies must have been con- 
taminated by the ejecta expelled during the AGB evolution of 
an early generation of intermediate mass stars (Cottrell & Da 
Costa 1981, D'Antona et al. 1983), or must have been formed 
by these ejecta (D'Antona et al. 2002). 

Recently, this scenario has been questioned by 
Denissenkov & Weiss (2001). Denissenkov & Herwig 
(2003), based on the computations of a 5M© model with 
metallicity Z = 0.0001, showed that a simultaneous depletion 
of oxygen and enhancement of sodium (as indicated by the 
anticorrelation observed) is made very difficult by sodium 
burning at high temperatures, those necessary to deplete 
oxygen; their analysis is independent of the amount of extra 
mixing assumed at the base of the convective envelope. 

We postpone to a forthcoming paper the discussion of the 
role which IMS may play in the pollution of the interstellar 
medium within GCs, because a more complete analysis, cover- 
ing the whole range of masses involved, and the most appropri- 
ate metallicities, rather than a single model, is required. 

In any case, the present work clearly suggests that among 
all the uncertainties included in the input physics adopted to 
calculate the evolutionary sequences the treatment of convec- 
tion has the greatest impact. AGB evolution is one of the few 
cases in stellar astrophysics where the results obtained are glob- 
ally dependent on the convective model adopted, in terms of 
time scale, luminosity, temperature and nucleosynthesis. 



Fig. 12. The variation of the surface sodium abundance with 
mass within the three models precedently discussed, and a fur- 
ther FST model where sodium production is favoured chosing 
an appropriate nuclear network. 

Let us summarize again what we can say about some of the 
important elements for the self-enrichment scenario. 

1 . Lithium: The constant lithium content of the ejecta, inde- 
pendent of the convective model, is in good agreement with 
the observed lithium abundances of NGC 6397 TO stars 
(Bonifacio et al. 2002), which all show approximately the 
same population II standard value. 

2. Helium: D'Antona et al. (2002) and D'Antona & Caloi 
(2004) have shown that the high helium content of the 
ejecta, confirmed by the present calculations, plays a role 
in the morphology of the horizontal branches of GCs. 

3. Sodium and Magnesium: From Fig.ll2lwe see that the MLT 
models can achieve sodium production quite efficiently, 
because of the TDU carrying fresh 22 Ne (later converted 
into sodium) to the surface; on the other hand, the FST 
model, living shorter because of a larger mass loss during 
the whole AGB evolution, ejects almost all the envelope 
mass well before the stage when TDU can operate signifi- 
cantly. 

In terms of sodium-oxygen anticorrelation, if the adopted 
mass loss rate holds, only the MLT21 model seems to pol- 
lute the interstellar medium in the "right" direction, be- 
cause apart from producing sodium efficiently it also de- 
pletes oxygen (see Table [2 and Fig. II It . though other 
problems arise, because in this case the polluted material 
would have magnetisium isotope ratios 25 Mg/ 24 Mg and 
26 Mg/ 24 Mg around ~ 3, which disagrees with the mea- 
sured magnesium isotopic ratios by Yong et al. (2003) for 
giant stars in NGC 6752. This problem is shared by the FST 
model. 

4. CNO: The C+N+O content of the ejecta is almost un- 
changed with respect to the initial value in the FST mod- 
els, which would be in agreement with the almost constant 
C+N+O value found for various GCs stars observed (M92: 
Pilachowski et al. (1988); NGC288 and NGC 362: Dickens 
et al. (1991); M3 and M13: Smith et al. (1996); M4: Ivans 
et al. (1999)). The MLT models lead to a strong increase of 
the C+N+O abundance, due to the effects of several TDUs. 

The self-enrichment hypothesis must be evaluated on the 
basis of a much deeper analysis, but it is clear that the uncer- 
tainties connected to convection can no longer be neglected in 
the computation of the ejecta of the various models. 

7. Conclusions 

We examine the effects of changing the treatment of convec- 
tion on the overall evolution of intermediate mass stars. We 
compare the stellar evolution models of initial mass M = 5M G 
calculated with the MLT treatment of the instability regions 
with those calculated with the FST model for turbulent con- 
vection. We show that while the evolutionary phases preceding 
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TPs are scarcely affected (with the exception of RGB and early 
AGB effective temperatures), the AGB evolution is heavily in- 
fluenced by the convective model adopted. 

A higher convection efficiency at the bottom of the outer 
convective zone of AGBs leads to a shorter distance between 
the edge of the CNO-burning shell and the base of external 
convection; this, in turn, leads to an extra- luminosity (which, 
in the first TPs, is entirely due to 3 He burning) which triggers a 
faster growth of the luminosity itself, hence of mass loss. 

From the physical point of view, the overall duration of 
the AGB evolution might be shorter by a factor of ~ 3 if the 
MLT treatment of convection is replaced by the FST model. 
Noticeable changes can be seen even if, still in the MLT frame- 
work, the free parameter a is changed. The maximum lumi- 
nosity achieved by the FST is ~ 70% larger, and the number 
of TPs (hence, of third dredge-up episodes) experienced by the 
FST model is considerably lower. 

As for the chemical content of the ejecta, the only robust 
conclusion which can be drawn is that these stars pollute the 
interstellar medium with material which is extremely Helium- 
rich (Fqecta = 0.31 in all cases) and with a lithium content 
which is only a factor of two lower than the standard popu- 
lation II value. 

The extent of oxygen depletion, of sodium enrichment and 
of the total C+N+O abundance of the ejecta are strongly de- 
pendent on the convective model adopted. 

Convection is one of the most important parameters in the 
AGB evolution. Any evaluation of the validity of the self- 
enrichment scenario for GCs must deal with the uncertainties 
connected to convection before deriving conclusions. 
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